Fractionalization via Z2 Gauge Fields at a Cold Atom Quantum Hall Transition 
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We study a single species of fermionic atoms in an "effective" magnetic field at total filling factor 
I'f — 1, interacting through a p-wave Feshbach resonance, and show that the system undergoes a 
quantum phase transition from a Vf = 1 fermionic integer quantum Hall state to = 1/4 bosonic 
fractional quantum Hall state as a function of detuning. The transition is in the (2 + 1)-D Ising 
universality class. We formulate a dual theory in terms of quasiparticles interacting with a Z2 gauge 
field, and show that charge fractionalization follows from this topological quantum phase transition. 
Experimental consequences and possible tests of our theoretical predictions are discussed. 



Quantum number fractionalization in condensed mat- 
ter systems is usually associated with topological order 
of the ground state. Although topological order inher- 
ent in fractionalization has been proposed in a number 
of systems including high- Tc materials [l|, Q and spin 
liquids [l, 0] , the only concrete example is in fractional 
quantum Hall (FQH) effect Q. Ultra cold atoms pro- 
vide an arena where topological order can be realized 
in quantum many body states by manipulating inter- 
actions. Interactions can be controlled by tuning the 
system through Feshbach resonances [B-Q. This allows 
one to investigate weak and strong paring regimes of 
fermionic superfluids with both s- and p-wave pairing in- 
teractions; the celebrated BCS-BEC crossover was clearly 
demonstrated for the former whereas the latter is ex- 
pected to yield a rich phase diagram with a variety of 
phase transitions as a function temperature and interac- 
tion strength Q . At the same time it has been suggested 
that FQH states can be realized in cold atoms by intro- 
ducing "effective" gauge fields through rotation of the 
atomic trap or atom- light interaction jlfl]. Abrikosov 
vortex lattices have been observed in these systems both 
through trap rotation 0] and generation of synthetic 
gauge fields [lO| . At sufficiently rapid rotational frequen- 
cies or subsequently low enough filling factor the single 
particle energies become highly degenerate and the sys- 
tem resembles an interacting electrons gas in a magnetic 
field; at this point the vortex lattice is expected to melt 
leading to various strongly correlated FQH states cor- 
responding to different filling factors. In fact FQH like 
correlations have already been reported in an ensemble 
of rotating traps 

The possibility of fermion pairing in an external mag- 
netic field [12I 1 allows one the unique opportunity to study 
Quantum Hall (QH) transitions driven by attractive in- 
teractions [13I ll^ . In this paper we consider fermions 
interacting through a p-wave Feshbach resonance in an 
external magnetic field at total filling factor Vf = 1. In 
the absence of interactions the fermions form an inte- 
ger QH (IQH) state, which is stable against weak par- 
ing interactions. In the strong paring (or "BEC") limit 
two fermionic atoms combine to form a bosonic molecule 
with twice the effective charge 2e*, and half the density 



of fermions; also the Landau level degeneracy is doubled 
due to the doubling of "charge", giving a bosonic fill- 
ing factor i/f, — 1/4. These molecules thus can form a 
FQH state of the Laughlin type for repulsive interac- 
tions. Although both states have the same Hall con- 
ductance a^y = {e*f/h = (l/4)[(2e*)2//i], they ex- 
hibit different topological orders and must be separated 
by a phase boundary. We use Chern-Simons-Landau- 
Ginzburg (CSLG) theory [l^ for QH states to show that 
this quantum phase transition (QPT) can be of second 
order in the (2 -|- 1)-D Ising universality class. To reveal 
the topological nature of the phases and phase transition, 
we perform a duahty transformation to show that the 
low-energy theory of the system is a Chern-Simons/Z2 
lattice gauge theory, coupled to a massive quasiparticle 
field. The transition occurs in the Z2 sector, which is also 
a confinement-deconfincmcnt transition for the quasipar- 
ticles. As a result quasiparticle charge fractionalization 
and corresponding change of statistics angle accompany 
this transition, as well as changes of other topological 
properties like ground state degeneracy on high genus 
Riemann surfaces [l6| . 

To study this QPT our point of departure is the CSLG 
theory, which relies on the off-diagonal long range order 
property of Laughlin type FQH states 17[ by mapping 
them to charged superfluids coupled to a Chern-Simons 
(CS) gauge field. In the weak paring limit we can view 
the fermionic IQH state in the following way: attach one 
quantum of CS statistical flux to each fermionic atom 
thereby transforming it into a composite boson; at the 
mean field level the flux attached to the fermion cancels 
the external magnetic field, and the composite bosons see 
zero net fiux. The fermionic IQH state then corresponds 
to the Bose condensed state of the composite bosons with 
charge-e* and can be regarded as an atomic superfiuid 
(ASF) albeit with a CS interaction. In the strong pair- 
ing limit two charge-e* bosons form a charge-2e* boson 
bound state; each charge-2e* boson carries two units of 
CS fiux quanta and the bosonic FQH state corresponds 
to the charge-2e* boson condensate and can be regarded 
as a molecular superfiuid (MSF) with a CS interaction. 
As we show below, this effective CSLG field theory also 
properly describes the QPT between these two phases. 
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Taking h = c = e* = 1, the imaginary-time action of 
the theory can be written as 
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Here ip^- are the bosonized atomic (a = 1) and the molec- 
ular (a = 2) fields, /io- is the respective chemical potential 
with IjL2 = 2^1 —S where 6 is the detuning; positive detun- 
ing favors atoms while negative detuning favors molecules 
with the Feshbach resonance corresponding to S = 0. 
m is the mass of the atom and g is the atom-molecular 
conversion matrix element, with the background density- 
density interactions implicit in " • • • " . The external mag- 
netic field is incorporated in the vector potential A, while 
a is the CS gauge field that attaches one (two) units of 
flux quanta to the atomic (molecular) fields respectively. 
The statistics angle 6 = tt encodes the flux attachment 
to the fermionic atoms transforming them into composite 
bosons. The CS statistical gauge field (restoring funda- 
mental units 00 = hc/e), 
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is chosen to cancel the external magnetic field V x ( (a) + 
A) = on average; thus all bosons experience zero net 
fiux. 

The action ([T]) defined above is invariant under the 
following local gauge transformations: 

^1 ^ ^ie^^(^'*); i2 ^ V2e2'^("^*); ^ ^ d^e{x,t). 

(3) 

The IQH state can be viewed as a charged Bose conden- 
sate of atoms anrf molecules: (i/'i) ^ and (V'2) 7^ 0; thus 
the local U{\) symmetry is "broken" (in an abuse of the 
term; see later for clarification) resulting in a mass to the 
gauge fluctuations Aa^ — a^+A^ through the Anderson- 
Higgs mechanism. The FQH state can be viewed as a 
charged Bose condensate of molecules (^^2) 7^ while the 
atoms have not as yet condensed: (V'l) — 0, again result- 
ing in a mass for Aa. Physically these reflect the fact the 
system is in an incompressible QH state throughout the 
phase diagram and charge fluctuations (which give rise 
to Aa) are suppressed. 

Despite the ordering of -02, in the FQH state a discrete 
Z2 degree of freedom remains unspecifled due to the tt- 
periodicity of 02 in ([3]). This is physically due to the fact 
that molecules are made of two atoms; thus specifying 
the phase of 02 only determines the phase of "01 modulo 
TT. Hence the FQH state only breaks the t/(l)/Z2 sym- 
metry spontaneously, while the IQH state in which the 



phase of tpi orders further breaks the remaining Z2 sym- 
metry. Based on symmetry grounds it is easy to antic- 
ipate that an Ising-like phase transition must intervene 
between the v\j = 1/4 FQH and Vf ^ 1 IQH states in 
which this hidden Z2 symmetry is spontaneously broken. 
In the following we show that this is indeed the case and 
derive the effective theory for the phase transition. 

To work out the effective theory describing the transi- 
tion we start in the vi, = 1/4 FQH region of the phase di- 
agram where as we have already mentioned the molecules 
Bose condense (02) ^ (assumed real and positive with- 
out losing generality) thereby assigning a Higgs mass to 
the gauge fluctuations: ((■02)^/4to)| Aap. This allows 
us to safely integrate out the massive CS gauge fluctua- 
tions Aa^. To make the Ising degree of freedom trans- 
parent we perform a transformation on the resulting ex- 
pression writing in terms of real scalar fields ipu and 0/ 
with 01 = tJjR + iipi, 



Ceff = 2iijRdriJl + Tr^|V0fl|2 ^ ^|V0^|2 
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where h — 5('02)- Since + 2h > fii — 2h, 0/ reaches 
criticality before 0^, the latter remains massive at the 
transition point and therefore can be safely integrated 
out of ^ to give us the effective theory for the transition. 
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which belongs in the (2 -|- 1)-D Ising universality class. 
This effective theory resembles the neutral atomic super- 
fiuid to molecular superfluid transition [l^ addressed by 
others in ultra cold atoms. The = 1/4 FQH to Vf — \ 
IQH phase transition is driven by the condensation of 0/ 
which happens at 8c —Ah; note that the presence of 
the molecular Bose condensate requires that fii > S/2. 

While CSLG theory correctly describes the phase 
transition using the language of conventional Ginzburg- 
Landau theory (in terms of broken symmetry and or- 
der parameters etc), it is important to note that there 
is no real symmetry change between the z/;, = 1/4 FQH 
state and the vf = 1 IQH state. This is because the 
order parameters of the CSLG theory, 0i and 02, are 
not gauge-invariant quantities. In fact this is a phase 
transition between two topological phases, with different 
topological quantum numbers. For example, the ground 
state degeneracy on the torus is 1 for the I'f = 1 IQH 
state whereas = 1/4 FQH state has degeneracy 4. To 
capture the topological nature of this QPT and correctly 
describe the change of topological properties of the sys- 
tem we must perform a duality transformation [li] on 
©. 

To formulate the dual theory we work with a discrete 
version of (IT|) on a square lattice, the lattice choice is 
arbitrary and we only require that the lattice version 
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correctly reproduce ([T]) in the continuum limit. The lat- 
tice hamiltonian can be expressed in the rotor formal- 
ism in terms of noncommuting [Oi^a-jfi'j.a] ~ i^ij number 
fii^c! and phase 9i^cr operators. The imaginary-time parti- 
tion function is obtained in the usual way by discretizing 
imaginary time into M — 1 steps of At — fi/M . Inserting 
the closure relations associated to the phase and number 
fields at each time step we arrive at, 
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where to- is the nearest neighbor hopping parameter for 
a = 1(2) atoms (molecules), the integers Ua^i repre- 
sent the atomic or molecular boson number and = 
1/2 + {pL„)/X and A > is the onsite repulsion. The lat- 
tice derivative is defined as A^O^.i = O^-i+i^ — 9a-i with 
6'o-;i a 27r periodic angular variable and i — {iT,ix,iy) 
representing sites on the spacetime lattice, a = x,y cor- 
responds to spatial directions and n = (r, a) denotes the 
spacetime directions. The lattice version of the interac- 
tion and CS term is. 



Sint = gAr ^cos(26'i;,;-6'2;i); Scs = ^^a^e^,.AA^ 
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To perform the duality transformation it is convenient to 
make the Villain approximation for the cosine terms in 
for example 



exp [t cos{AaOi)] 



exp[ty(A„0, -27r;„;,)'], (9) 



with tv a renormalized constant. Since ^ preserves the 
periodicity of 9 this change in functional form should not 
alter the physics. Hereafter we denote the modified par- 
tition function as Zy and drop the unimportant normal- 
ization constants without comment. Following standard 
manipulations [191] we introduce Hubbard-Stratonovich 
fields J^.j, J^.j, r]i which allows us to integrate out 9^.1 on 
all sites and intermediate times to arrive at the current 
representation of the Villain approximation of (jH]), 
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where {Lf,,,} = {L^.i, Ll.i,r],} with L"^.^ = J^.^ - n„5^,r- 
J^.j = [n^-i, J^j, J^.j) is the integer boson 3-current and 



we have chosen At such that taAr = l/(AaAT). We 

transform to new link variables r,„ = j},, + 2J? and 

ifi ifi 

Kifj_ — J}^, modifying the constraints to A^r^;i = 
and Afj_Kfj_-i — 2?/;. The first constraint physically corre- 
sponds to the fact that the total atomic and molecular 
current F^^ is conserved whereas the second constraint 
encodes the local atom-molecule interconversion events. 
The second constraint A^if^ji = 2r]i can be used to sum 
over rji. We can further sum over the link variables K^-i. 
This summation is constrained as K^-i and F^;i have the 
same parity. To keep track of this parity, we divide F^.^ 
into an even part and an "Ising" part which are indepen- 

+ Here F;^., is even 



dently conserved: F^. 



while F^.j can only take values —1,0,1. This is always 
possible due to the total current conservation; links with 
odd Fp-i must form loops, which can in turn be identified 
as loops of F^.j ^ 0. Therefore in the summation if^.^ 
has the same parity as that of F^.^. After this summa- 
tion, the partition function in terms of F"^ and F^.j takes 
the form, 



Zv= Va^ J2 exp ( XI [ 

•' r= ^ i.u 
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- 27rjF^;i(5a^;i] + ^ ^ a^^^e^^^yxA^ax H j , (11) 

with the constraints A^F^.^ — A^F^.^ — 0. They can be 
solved as usual by defining integer link variables a and 
/i^t;a on a dual lattice, 

^iti — C/ii/AAixf/Aia, — ^iii^xAi,hx;a, (12) 

where a denotes the sites of the dual cubic lattice and 
.g^i;a is an even-valued integer field and we define h)j,-a as 
an integer valued field. To proceed we remove the even- 
integer constraints on Qfj^-a by performing a Poisson re- 
summation through half-integer link variables m^-a and 
exploit the gauge invariance Trg^^a 7i'5/i;a ~ A^ipa to 
introduce tpa which we later identify as the phase of the 
molecular vortex field. This choice is slightly different 
from convention [l^, and reflects the evenness of g^-a- 
Adding a vortex fugacity term l/yJ2^i-a we can 

integrate out m^-a at the expense of another Possion re- 
summation even-integer dual link variable n^;^ to obtain 

Zv = j VipaVa^Vg^ ^ exp ^ ^ [^a^e^^AA^c 

-I- 2Trief^^xA^{gx-a + hx-a)Sai_,.i - /3 cos(7re^^AA,,/iA;a)] 
+ yJ2{2Af,ipa - 2ngf,,a + 2Tm^,,af + • • • ) (13) 



where (3 — 1/ {AArt'i ) , from here onwards the Maxwell 
and other higher order terms will be implicit in 
Redefining the U{1) gauge field 6^:a = 5^:0 -I- h^j^-a in the 
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resulting expression and integrating out the CS gauge 
field which is coupled to 6^;a (whose curl is the total con- 
served current), one arrives at the dual version of ([5]): 



^ J VipaVb^ ^ exp ^ cr^-a cos(A^(^a - 7r6^.a 



□ 



(14) 



where we have introduced the Ising gauge field 
f/iia = cos(7r/i^;(j). In (jl4p the angular variable (/Jq 
corresponds to the phase of the molecular vortex field. 
Due to the presence of the (dual) CS term such a 
vortex carries charge e*/2 = 2e*/4 with 7r/4 statistics. 
It corresponds to the quasiparticle excitation of the 
i^b = 1/4 bosonic FQH state. 

In d -|- 1 > 3 Z2 gauge theory exhibits a confinement- 
deconfinement phase transition as a function of (3 [20| . 
The difference between the confining and deconfining 
phases relies on the notion of topological order which 
has important implications on the nature of excitation 
of the system. Pertinent to our case is the charge 
fractionalization that occurs in the deconfined phase: 
e*/2 quasiparticles (molecular vortices) are free to prop- 
agate, this corresponds to the molecular FQH phase. In 
the confined phase the energy grows linearly with the 
separation of two e*/2 quasiparticles: they thus bind 
together and form charge-e"^ with 4 x 7r/4 or fermionic 
statistics, these are nothing but the original fermionic 
atoms. Thus the confined phase corresponds to the 
atomic IQH phase. A more physical way understand 
this transition is in terms of condensation of visons, 
which are excitation of the Z2 flux through a plaquette 
(Tin''' — 0- the confined phase visons are 

condensed and particles carrying Z2 charge cannot 
propagate, whereas in the deconfined phase where the 
visons have not as yet condensed these topological 
excitations cost finite energy, and the theory is in the 
deconfined phase. Furthermore (|14p correctly captures 
the topological degeneracy of the phases on high genus 
surfaces. Using torus as an example, this degeneracy 
is 4 in the deconfined phase as in the Vf — 1/4 FQH 
state, and 1 in the confined phase as in the Vf = 1 
IQH stated, Thus the dual theory (HH) captures all 
the topological properties of the phases, and properly 
describes the phase transition from i^b = 1 IQH state to 
1/^ = 1/4 FQH state. 

For d — 2 this confinenient-deconfinement phase 
transition is in the 3-D Ising universality class [20| . 
The critical exponent of the 3-D Ising model are well 
known and standard scaling arguments 22| predict that 
a gap vanishes as one approaches the critical point: 
A ~ |(5 — dc\^'^ , where i/ !=s 0.63 is 3D Ising correlation 
length exponent and the dynamical exponent z = 1 due 
to (emergent) Lorentz invariance. Physically this gap 
corresponds to the long-wavelength magneto-phonon 



gap in the Vb = 1/4 FQH state which can be detected 
by probing the excitation spectrum through simulated 
Bragg spectroscopy [23| . 

The stability and observability of the FQH phase 
and the critical point depend on the repulsive nature of 
the background molecule-molecule interaction assumed 
here, and in particular, the stability of the atoms 
and molecules themselves. While in 3D their decay 
rate is comparable to the interaction scale near the 
p-wave resonance, the situation becomes much more 
favorable in quasi-2D systems The presence of 

an external "magnetic field" may further help by pre- 
venting multiple particles from coming close together. 
We also note that the transition occurs close to, but 
not precisely at the resonance; the distance between 
criticality and resonance depends on density and other 
parameters. This allows for additional room in searching 
for a parameter window where the critical point is stable. 
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